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Abstract 

We consider an infinite system of non overlapping globules undergoing Brownian mo- 
tions in R 3 . The term globules means that the objects we are dealing with are spherical, 
but with a radius which is random and time-dependent. The dynamics is modelized by 
an infinite-dimensional Stochastic Differential Equation with local time. Existence and 
uniqueness of a strong solution is proven for such an equation with fixed deterministic ini- 
tial condition. We also find a class of reversible measures. 

AMS Classifications: 60H10, 60J55, 60J60. 

KEY- WORDS: Infinite-dimensional Stochastic Differential Equation, hard core potential, 
oblique reflection, reversible measure, local time. 
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1 Introduction 



The aim of this paper is to construct a random dynamics performed by an infinite system of 
globules, where a globule is a sphere in M 3 with variable radius. The centers of the globules 
undergo independent Brownian motions, while their radii perform Brownian oscillations 
between a minimum and a maximum value. Since the scale of these oscillations can be 
different than those of the centers, we introduce a coefficient a which reflects the elasticity of 
the surface of each globule. The globules can not overlap and when the distance between two 
globules becomes 0, they repel each other immediately; that means they interact through 
a hard core potential. 

A reversible system of infinitely many Brownian hard spheres (called balls) was first 
introduced and analyzed by H. Tanemura [UJ. Then, some natural generalizations were 
studied for different types of additional smooth interactions between the balls: for a gradient 
type interaction with finite range in [I], [5] and for an interaction with infinite range in [7]. 
The specificity one has to deal with in a hard core situation - hard balls can not overlap - 
comes from the additional infinite-dimensional local time term in the SDE. Notice that in 
all these works the spheres have a fixed positive radius. 

The originality and new difficulty of the present model - which can find relevant ap- 
plications in cell dynamics like molecular motors - lies in the random oscillations of the 
radius of each sphere. We propose here a pathwise approach for the construction of this 
infinite-dimensional dynamics, by building a sequence of finite-dimensional approximating 
processes. But already for finitely many globules, the existence of such dynamics is not a 
simple question. Indeed, one of the authors constructed recently in [3] a finite system of 
mutually repelling Brownian globules. Nevertheless, we need here a non trivial generaliza- 
tion of these results: since the scale a of the radii oscillations is different than the scale of 
the center oscillations, the direction of the reflection after a collision between two globules 
is no more normal as in [3]. It is an oblique reflection of Brownian motions on a complex 
non smooth domain, whose existence problem we solve in Proposition 13.11 

In Section [2] we present the model and its dynamics described by the stochastic differ- 
ential equation (£) and we state the results. In Section [31 we show the convergence of the 
approximations and analyze the limit process. Last, we remark that some kind of hard 
core Poisson measure is reversible for this dynamics. 

2 The infinite model of mutually repelling globules 
with Brownian radii 

A globule is a sphere in M. d with a variable radius. For d = 2 the globules modelize for 
example the motion of discs on a flat surface or balls floating on a liquid. In this paper, we 
fix d = 3 which corresponds to the natural physical case of bubbles in the Euclidean space. 
Our techniques and results obviously extend to any dimension d larger than 1. 

A globule is characterized by a pair (x, i) G K 3 x R. x is the position of the center of 
the globule and x is its radius. 

We are dealing here with infinitely many indistinguishable globules, thus the state space 
of the system is included in A4, the set of point measures on M 3 x R. A configuration 
of globules is a locally finite point measure x = YlieJ^fafo) on ^ 3 x ^> wriere ( x ii%i) 
characterizes the i-th globule and J C N. For simplicity sake, we will identify any such 
point measure x with its support {(xj, xi), i G J} C M 3 x R. 
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The globules we deal with in this paper can not overlap and their radii are bounded 
from below (resp. from above) by a constant r_ > (resp. r + > r_). Thus the exact 
configuration space A g of all allowed configurations of globules is the following set: 

. f x = {(xi, Xi),i E J} for some JcM, 1 
9 with Xi E IR 3 , Xi E [r_; r + ] and |xj — > Xj + for i ^ j J 

Let us notice that this model can not be reduced to a hard core model in IR 3 x R. In- 
deed, a hard core condition between globules would mean that there exists p such that 
Vz ^ j, \(xi,Xi) — (xj,Xj)\ > p, which is equivalent to \xi — Xj\ 2 + \xi — Xj\ 2 > p 2 . This 
last inequality is clearly not comparable with the condition | 

We will use the notations: 

• B(x,p) is the closed ball centered in x E M 3 with radius p and by extension, for any 
subset A in IR 3 , we define the p- neighborhood of A by 

B(A,p) := {y E M 3 such that d(y, A) < p}. 

where d(y, A) denotes the Euclidean distance between y and A. 

• The symbol \v\ denotes the Euclidean norm of the vector v. 
We also denote the volume of a subset A of IR 3 by | A | . 

• For A x / a Borel subset of M 3 x IR, Naxi is the counting variable on Ai : 

N AxI (x) = (J{i E N : Xi E A and x { E I}. 



• For A a Borel subset of IR 3 x IR, B\ is the u-algebra on Ai generated by the sets 
{N A , = n}, n E N, A' C A, A' bounded. 

• We write xa = x fl A for the restriction of the configuration x to A C K 3 x 1, 
and xy for the concatenation of configurations x and y. 

• 7i (resp. 7Ta) is the Poisson process on I 3 x R (resp. on A) with intensity measure 
the Lebesgue measure dy (resp. dy\\). 

We define the set Il 9 of hard globule Poisson processes via a local density function : 

Definition 2.1 A Probability measure p on Ai is a hard globule Poisson process if and 
only if, for each compact subset A C IR 3 x IR, 

li(dx\B A °)(y) = 7^!{x A y A ceAj ^A(rfx) for p-a.e. y. 
where the so-called partition function Z A,y is the renormalizing constant : 

Z A ' y = e~ |A| [ 1 + ^2 ~T / ^AYAoeAg} dx x dx x • • • dx n dx n 



3 



At least one hard globule Poisson process exists (generalization of the existence results for 
hard core Gibbs measure in [I]). It is conjectured that for r + small enough it is unique, 
while for r_ large enough phase transition occurs : U g should contain several measures (see 

e.g. [H]). 

In order to modelize the random motion of repelling globules with oscillating radii, let 
us consider a probability space (f2, J 7 , P) endowed with a complete filtration {J-t}t>o and 
two sequences of J^-Brownian motions: (Wi(t),t > 0) ie ^ which are independent M 3 -valued 
Brownian motions and (Wi(t),t > 0)^ which are independent R- valued Brownian mo- 
tions, independent from the WiS too. 

We fix a parameter a > which measures the scale of the radii oscillations, that is the 
elasticity of the surface of each globule. 

We consider the following system of stochastic differential equations with (oblique) reflec- 
tion : 



For i e N,t e [0,1], 

Xi {t) = Xm + Wi {t) + W' 

7TZ Jo XAs) + Xj(s) 



jevt 



X i (t)=X i (0) + aW i (t) -o- 2 J2 L v (*) - L i+ (*) + L i- (*) 

where the local times L^L^L^ satisfy L^t) = / l {lXi{s y Xjis)l= x i{s ) + x j{ s)} dL 

t t 
U+it) = / l {Xi(s)=r+} dL l+ {s) and L;_(t) = / l {x . (s)=r _ } dL^s). 



As usual, the collision local times are non-decreasing M + — valued continuous processes 
with bounded variations and satisfy Ljj = Lji and La = 0. The starting configuration 
X(0) = {(Xi(0),Xi(Q)),i e N} is a point in A g . 

A solution of the system (£ g ) is a family ( Xi (t), Xi (t), Lij(t), L i+ (t), Lj_(t), < t < 
G N) of processes satisfying (S g ). 

Let us interpret the different terms of (S g ) : 

• when two globules collide (\ Xi (t) — Xj(t)\ = Xi(t) + Xj(t)), they are deflated (Xi(t) 
decreases by dLij(t)) and move away from each other (Xi(t) is submitted to the 
repulsive force 

• when the radius of a globule reaches the maximal value (Xi(t) = r + ), it is deflated 
(Xi(t) decreases by dLf(t)); 

• when the radius of a globule reaches the minimal value (Xi(t) = r_), it is inflated 
(Xi(t) increases by dL~{t)). 

Theorem 2.2 The stochastic equation (S g ) admits a unique solution with values in A g for 
any deterministic initial configuration which belongs to a full measure subset A g in A g . 

Proposition 2.3 If the initial distribution \x is a hard globule Poisson process and if 
[i(A g ) = 1, then the solution of (£ g ) is time-reversible, that is its law is invariant with 
respect to the time reversal. 

The next section is devoted to the proofs of these results. 
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3 The infinite-dimensional process, constructed by ap- 
proximation 

3.1 The approximating processes 

In this whole subsection, £ G N* and y G A g are fixed. Using a classical penalization method 
with external configuration y, we construct an approximating process which essentially 
stays in the ball B(0,£) . This is done by introducing in the dynamics (£ g ) an additional 
drift which vanishes in a subset of B(0,£) and is strongly repulsive outside of B(0,£); we 
take as drift the gradient of the C^-function (twice differentiable function with bounded 
derivatives) which is defined on R 3 x R by: 



ip e ' y (x,x) = ipi(\x\) +ip 2 (x) + ^ 3 

r-\vj\>e 



\ x - v_i 



with ipi, ip2 and ^3 non-negative C°°— functions vanishing respectively on ] — 00, 
and [1, +00 [, and increasing rapidly on their supports : 
V>i(s) = 2s for s > £ + e~ £ 

■02 (s) = £s for s > r + + e~ e and ij)<z{s) = £{r + + r_ — s) for s < r_ — e~ l 
ip3(s) = £ for s < 1 — e~ e . 
The function ip e,y satisfies 



r_, r 



■?//' y (£,£) = a: G B(0,£) and (x, x)yB(o,. 



eA g . 



By a slight abuse of notation, y B ( 0> f) C denotes the restricted configuration ^-i{j:\ yj \>£}{(yj, Vj)}- 
Remark that the functions ip e,y are so repulsive for large £ that they satisfy 

+00 „ 

sup Y / V ,y (*.*)>o exp(-ip e ' y (x,x)) dxdx < +00. (1) 



Let us now define the finite-dimensional dynamics : 
f Vi£{l,...,ra}, VtG[0,l], 



XAs) 







Xi(t) = Xi(0) + Wi{t) - - / V^' y (^i(3),Xi(a))da + J] 

2 /-t « 
= 1,(0) + aWi(t) -— V^iX^Xi^ds - a 2 Y,Lij(t) - L i+ (t) + U_(t) 

Jo .7=1 



where the local times satisfy Lij(t) 



L i+ (t) 



1 



{Xi(s)=r + } 



dL i+ (s) and Li-(t) 



1 



{Xi(s)=r-} 



dLi-(s). 



(£„) is a n-dimensional reflected stochastic differential equation. 

If a = 1, that is if the radii oscillations and the center oscillations are on the same scale, 
the above Skorohod equation contains a normal reflection on the boundary of the set of 
allowed configurations of n globules. The problem of existence and reversibility of this type 
of dynamics was recently solved by one of the authors in [3]. 

For o ^ 1 a new difficulty occurs. The physically evident reflection on the boundary 
of the domain of allowed configurations is now oblique since the radii oscillations have 
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another time scale as the center oscillations: it corresponds to a normal reflection but in 
an anisotropic configuration space, where the radii coordinates are rescaled by 1/a. The 
existence of solution for general SDEs with oblique reflections on nonsmooth domains is 
a hard problem which is solved in the literature only in some particular cases, like for 
intersection of smooth bounded domains or for polyhedral domains (see e.g. [2] and [12]). 
Since our model is not covered by these works, we present in the following proposition a 
suitable existence result. 

Proposition 3.1 Assume that <3> is an ^-valued C^-function defined on (M 3 x M) n . There 
exists a unique strong solution to the Skorohod problem 

{ Vie{l,...,n}, Vte [0,1], 

X<(t) = X,(0) + Wi{t) -\f V x MMs))ds + £ f ^^(s)dL lJ (s) 

z Jo „■_-, Jo Xi + 



(*») 4 



j=1 Jo Xi + Xj 

2 rt n 

Xi{t) = 1,(0) + &Wi{t) - °— / w £ MMs))ds -o- 2 J2 L a(t) - L i+ {t) + LUt) 

t j=1 

where the local times satisfy L y (t) = / ^x^-x^s^x^+x^s)} dL v( s ) » 

t ft 

L i+ {t) = / l {Xi{s)=r+} dL i+ (s) and L*_(t) = / l { x i(s)=r _ } dL^(s). 
Jo Jo 

For any initial condition in A g the solution is an A g -valued process. 

Moreover, the solution with initial distribution e - *^!^ (x)dx is time-reversible. 

The proof of this proposition is postponed to the end of this section. A key idea is the 
transformation (see fl2J)) of the initial Skorohod problem into a simpler one, by stretching 
the radii coordinates by the factor 1/a and thus transforming the oblique reflection in a 
normal one on a modified domain. Let us underline that the oblique reflection we consider 
is the unique one for which the existence of a reversible dynamics is ensured. 

Applying Proposition 13.11 with the potential $(x) = {xi, Xi), we obtain the 

existence of a solution to (S^ ,y ). 

When the initial condition is the deterministic configuration x, this solution is denoted 
by X £,y,n (x, •). In particular, the *A 9 -valued finite-dimensional process X^ y with initial 
configuration x = yB{o,e) and external configuration yB(o,e) c evolving under the random 
dynamics (S^' y ) is : 

X^(-) := X^ n {y B{0/) , ■) with n = ${i G N : Vi G B(0, £)}. 
The associated local times are denoted by Lf j, L^, , i, j G N. 

If the initial condition of the system (£^ ,y ) is random with distribution given by the 
finite measure: 

n 

u^ y (dx) := exp(— ip e ' y (xj, xj)) I^ 9 (x) dx\dx\ ■ ■ ■ dx n dx n . 

i=l 

then the solution of (£^' y ) is reversible. Its law is denoted by ' . 
Consider now the following Poisson mixture in n of the Q^' y 's: 

-\B(0,£)\ +2^ l 



n £ ' y = - V — o^ y 

v Zt,y n \ v « ' 

n=0 
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where Z i,y is the renormalizing constant which ensures that Q e,y is a Probability measure. 
As a mixture of ^-supported time-reversible measures, Q e ' y is time-reversible with support 
included in A g . Its projection at time is the Probability measure 



//' y (dx) 



1-8(0,- 



n=0 



which represents the law of a Poissonian number of globules essentially concentrated in 
B(0,£). 

We will construct the infinite-dimensional globule process as limit in i of X^ y ; unfortu- 
nately, X £,y is not time- reversible. This is why we had to introduce Q e,y , whose reversibility 
plays a crucial role in the study of the set of nice paths defined in the next section. More- 
over, we will prove that the law of X £,y and Q E ' y are asymptotically close. 

To complete this section, let us prove Proposition 13.11 

Proof We first introduce an anisotropic linear transformation er _1 on the space of globule 
configurations by 



Vi, xl = Xi and x° 



1 



-Xi 



a 



We also transform the process, the potential and the local times as follows : 



a 



■Li. 



a 



(2) 



X° = tr- 1 X, $"(■) = $(<t-), L% = V2 + 2^L lv L^ 
Moreover, the set of allowed configurations becomes 

Ag := {x : <rx G A g }. 

X and its associated local times are solution of the system (8 *) if and only if X CT and its 
associated transformed local times are solution of the following system: 



x?{t) = x° (o) + Wi(t) 



V^(X a (s))ds 



x°-x° 



[s)dLUs) 



a 



V2 + 2a 



V2 + 2a 2 a(X[ + X?) 

n 

"£LUt)-L° + (t) + LUt) 



where the local times satisfy L^(t) = L^it) = J i^ x? 



j"=i 



L° + (t) 



I 



- } dL° + {s) 



and L1_{t) 



1 



- } dLUs). 



Furthermore, the reversibility of X is equivalent to the reversibility of X ff ; with other words, 
the solution of with initial distribution ^e _ *^l^ g (x)dx is reversible if and only if 
the solution of with initial distribution -^e~^ a ^ y 'l^<r(y)dy is reversible. 

The new system of globules (£* ,<J ) has now the form of a Skorohod problem with normal 
reflection. Thus it has a unique solution under the assumptions of Theorem 3.3 (and 
Corollary 3.6) in [3], that is if the domain Ag on which the equation is reflected satisfies 
the geometrical regularity properties listed in [5] Proposition 3.4. The rest of the proof 
consists in showing these four properties (see [5] for the relevant definitions). 
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The domain A° does not have a smooth boundary but it is the intersection of smooth 
domains in the following way: 

a ° 9 = ( n p v) n ( n n ( n 

l<i<j<n l<i<n l<i<n 

where Vij = {x, \xi — Xj\ > a{xi + Xj)}, T> i+ = jx, Xi < — j and Dj_ = |x, > — |. 

(i) At each point x 0/ £/ie boundary of the smooth set T>^ (resp. T> i+ , T>j_), there exists a 
unique unit normal vector n^x) (resp. rij + , rij_,). 

Each T>ij is a smooth set with unit inward normal vector at point x e c^Dij equal to 
rijJx) = = where w« = —7- ^-r- = — if; = Wj = —a and wu — wu — 0, k % 

Each T>i + (respectively T>i ) is a half-space of (M 3 xK)" with a constant unit inward nor- 
mal vector n i+ = (0, ■ ■ • , 0, —1, 0, ■ • • , 0) ((2i — l) th coordinate equal to —1) (resp. rij_ = 
-n i+ ). 




Figure 1: A configuration of 5 globules in &D\- fl dV 2 + H <9£>3 ; 4 R and the different 
directions of impulsion ni__, n 2+ or n 3;4 (x) to go back into the interior of A°. 

(ii) Each set has the Uniform Exterior Sphere property on A a g : 

3a tj > 0, Vx e Ag D 92? y 5(x - n (; u, ; (x). ay) n % = 0. 



Each x e AZ fl <9Py satisfies — x 



a(xi + Xj) > 2r_ 



For x es = x — r_ \Jl + 2cr 2 rijj (x) and for any z : 

I (xf + Zj) - (o^ s + Zj) I - cx(ij s + Zi + xj s + %) 



< 1^1 + |zj| - a{zi + %) + 



.1 j 



Xi Xj 

Xj — 2r_ 



< V2 + 2a 2 \z\ - 2r_(l + a 2 ) 



cr(a:; + 



— cr(xj + Xj) — 2r_cr 



This is negative as soon as |z| < r_V2 + 2cr 2 . Consequently, the property (m^ holds with 
ay = r_A/2 + 2<t 2 . See figure 2. 
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Figure 2: A configuration x G dT),^ and the corresponding center x es of the Uniform Exte- 
rior Sphere : x es = x — a^-n^x) = x — r_ \Jl + 2<r 2 rijj (x) . Left, a simplified representation 
in IR 2 . Right, a representation as a pair of colliding globules. Remark that x es (jL ^4^. 



fmj i?ac/i set T>ij has the Uniform Normal Cone property on A°: 

Bflij G [0, 1[ and 5ij > such that, for each x G A a g fl dT>ij, there is a unit vector 1^ 
satisfying 



y G A* n 5Z? y n 5(x, <y =► n«(y).15 > ^/l - 

For x,y G ^ n <9£>^ one has n^(y).n ij (x) = — ^ ( I % I ^'i + ff2 ) ' The 
inequality 



■''-■'•'.^ * >* V ^ 1 ^ 1 >l-2y^ |x ~ y| >l-^|x-y| 
- ^"1 - Vj\ l + V^j^L \Xi-Xj\ 



r. 



implies that ny(y).iijj(x) > 1 ^ ; _ 9 J x — y|. Hence, for 1^ = riy(x), (Hi) is satisfied 



V2_ 
r_(l + o 2 ) 

for each /3 i<7 - G [0, 1[ as soon as 5^ < — ^J* (1 — \J^~ Pfj)- 
(iv) Compatibility between the boundaries: 

( xe9Hy=^v(x).n ij (x) >/?o|v(x)| 
3/3 > y / 2max~~^~, Vx G cMg, 3 v(x) stzc/i £/iat < x G 9Dj+ =>- v(x).n i+ (x) > /3 |v(x) 

[ x G 9Z?i_ v(x).n 4 _(x) > /9 |v(x) 

Let us define the following cluster 

Cj(x) = {z} U {A; m , 3fci, ...,k m such that x G &D ikl n <9£> fclfe2 D • ■ ■ H <9P fcm _ lfcm } 
and define the vector v(x) = (fi(x), £>i(x), ■ • • , f n (x), #n( x )) by 

Vi G {1, • ■ ■ , n} t,,(x) = - ^ Jj^ x fc and *(x) = ^—^-Z^ r _. 
See figure 3. 

Since |i>i(x)| < (n — l)(2r + ) and |#j(x)| < 2 {&vi) ' then l v ( x )l — 2r + n 3 / 2 . Moreover : 
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Figure 3: A cluster with 3 globules and the associated impulsion v(x) constructed to push 
x G &D back into the interior of the set of allowed globule configurations. 



• if x G dVij, then C$(x) = Cj(x) and 

\/2 + 2a 2 v x ).n lj x = a x* + x 3 ) - ar_^- \ J > 2r_ - — — > r_ 

(r + — r_)(cr V 1) a V 1 

• if x G <9£> i+ , i.e. aXi = r+, then v(x).n i+ (x) = 

• if x G &Di-, i.e. ax; = r_, then v(x).n»_(x) = 2 (gvi) • 

r_ 

So, with /3 = -; — t , , , (iw) is satisfied. ■ 

4r + (cr V ljri^/ 2 

3.2 A full set of nice paths 

From now on, the techniques we use to study the globule model present some similarities 
with the methods developed for the model of hard balls treated in [5]. So, in the rest of 
the paper, we will only detail the proofs which contain new technical difficulties. 

We first bound from below the probability of globule paths which do not move too fast 
under the (£^ ,y )-dynamics. 

For every e > and 5 g]0, 1], let Af(5,e) denote the paths for which all globules have a 
(5-modulus of continuity w smaller than e, i.e. 

Af(5,e) = \xeC([0,l],A g ): Vf, w((Xi, 5) < e} , 
where the 5-modulus of continuity of a globule path (X,X) on [0, 1] is defined as 

w((X,X),6):= sup J\X(t)-X(s)\> + (X(t)-X( S )y. (3) 

0<s,t<l 
|t-s|<<5 



Proposition 3.2 There exists c > and c\ > such that the following lower bound holds : 
V£>0,V<5 G]0, l],WeN*, 

£ 3 ( e 2 



inf Q^ y {Af{5,e)) > 1 - c x — exp -c s 
yeAg \ 
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Proof of Proposition 13.21 

By construction, the processes 

Wi(t) = X^' n {t) - Xf y ' n (0) 

2 ^Jo J 5 ^' (s) +X/ y ' (s) 

and^(t) = i(^ B (t)-^(0)) 

(T /"* ™ 1 1 

+ - / V^' y (^ y ' n ( 5 ),^( S ))rf S + a^^(t) + -L, + (t)--L,_(t) 
2 J a a 

are 3-dimensional (resp. 1-dimensional) Brownian motions starting from 0. 

When the initial distribution is the law Q 1 ^ of X £ ' y,?1 is reversible, and the backward 

processes 

W t {t) = X*> y > n {l-t)-Xf' y ' n {l) 

1 r 1 n r 1 x e ' y ' n (s) - x e ' y ' n (s) 

2 Ji-t Ji-t X i y ' (s) + X/ y < (s) 

and^(t) = -(l^ re (l-0-^' y ' n (l)) + f T V*il&{Xt*> n {8),%*> n (8))da 
° 2 Ji-t 

n 

+a - L^l ~ t)) - -(Ml) - L l+ (1 - t)) + -(Ml) - - *)) 

3=1 ^ ^ 

are Brownian motions too. 

As in [9], the above equations provide the identities 

X e,y,n (t) _ x i,y,n {Q) = 1 ^ + ^ (1 _ t) _ 
Xf'^t) ~Xf' y ' n (0) = ^{WM + Wiil -t) -Wi(l)). 

Therefore, the control of the modulus of continuity of a globule path (Xf ,y,n , Xf ,y,n ) reduces 
to the estimate of the modulus of continuity of Brownian paths, as follows: 

Q £ n y (meY) 

< 2 n P (w((Wi,<7Wi),<?) > e) z/' y ((K 3 x R) n ) 

< 2 n P (w((Wi,c7Wi),5) > e) [ exp(-^ ,y (xi,zi)) ^c/xi ^((M 3 x M) n_1 ) 

Now, we can use the following estimate obtained as corollary of Doob's inequality (for a 
proof in the one-dimensional case, see the Appendix of [5]) : 

Lemma 3.3 Let us consider two independent Brownian motions B E M 3 and B E R. 
There exist two constants c > and C2 > ( depending only on a) such that for every e > 
and every S e]0, 1] 

c 2 / £ 2 



P(w((B,aB),5) >e) < ^expl- Cj 
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This leads by summation in n to : 

-\B(0 i)\ +°° i 
n=l 



— lB(0i)\ +°° i o 2 f 

^ e -^r{Hl TY/n*im 3 x Rf- 1 )) ^exp(-4) / exp(-^(x 1 ,x 1 )) dxid*i 

< — r^exp(— c^-) / exp(— ip e ' y (xi, xi)) dx\dx\. 



! x [r_,r + ] 

Recalling that ^r ,y only vanishes into the ball B(0,£), we get 



exp(—ifj e,y (xi,xi)) dxidxi 

J B(0,e)x[r-,r+] J 

The first term of the right side is smaller than £ 3 (r + — r_)\B(0, 1)| and, thanks to (JTJ), the 
last term is uniformly bounded in i and y. This completes the proof of Proposition 13.21 ■ 

In order to control the convergence of the finite-dimensional systems, we have to estimate 
how many globules collide with a fixed globule % during a short time interval. If the paths 
have a small oscillation, this set will be finite because globule % can not reach globules which 
are too far away. But we also have to avoid the bump to propagate along a large chain 
of neighboring globules. We first define patterns called chains of globules, and then prove 
that they are rare enough, in the sense that their probability decreases exponentially fast 
as a function of the length of the chain. 

Definition 3.4 Let e > 0. The set of configurations containing an e-chain of M globules 
is defined by: 



(^) \^ ^ *^~9i ^^1? ' ' ' i^M distinct , |^^2 *^%\ | ^ ^ ^ i 5 I^^m X ^M— 1 1 ^ ^ 



I'M <~ X iM-l~ 



We now define a set of paths which are smooth in the sense that, at regular time intervals, 
there is no chain of globules: for 5 G 1/N*, M G N*, e > 

tf(6,M,e) := |xgC([0,1],^ 9 ) : Wk e {0, . . . , ^ — 1}, X(5k) £& M (e 

Note that this set decreases function of e. 

We now prove a lower bound for the Q^ y -Probability of N(S, M, e). 

Proposition 3.5 For any M G N* ; there exists C3 > such that, for any 5 G 1/N* and 
< e < 1 : 

inf Q l * (Sf(S,M,£)) > \-%e £ M -\ 
y&A g \ J 
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Proof of Proposition 13.51 

Let us estimate the z/' y - and the //^-Probability that a chain exists. For n > M : 



M M 



/ n' y ( CI)A/ ( £: )) — 7 7777 / T lxi+Xi- 1 <\x i -Xi- 1 \<Xi+Xi- 1 +e ]Jlr_<Xi<r + 

(n-M)! j (r3xr)™ ~_2 7=7 

AI n 

exp(-^^ £ ' y (xi,Xi))exp(- ^ ^ y (x u Xi)) 

i=l i=M+l 

lr_<Si<r + Izi+Xj^ixi-Xji dx\dxi ■ ■ • dx n dx n 

i=M+l M+l<i,j<n 

f 

R 3 xR" 



5(0, xi + x 2 + e) \ 5(0, xi + x 2 ) | ■ ■ • (5(0, + xjw + e) \ 5(0, x M -i + x M ) 
J I r _<^< r+ exp(-^ £,y (xi, xi)) dxi cZ£i • ■ ■ cZx M 



i=l 



< (n ^ M) , ^M((^ 3 xR)- A/ ) |5(0,2r + + e)\5(0,2r + )| A/ - 1 



exp(— t/j ' y (xi, Xi)) dxidoti 

•X [r_,r+] 



'.y CCTB>3 v TO^-M\ „_*3,Jlf-l 



for a certain constant c 3 > 0. Therefore 

n=M ^ >' 

and the stationarity of Q e ' y implies 



i — n \ / 



k=0 



To prove the convergence of the approximations, we have to connect in a right way the 
different parameters 5, e, M, £, in order to introduce a set of nice paths Q y C Q on which 
the convergence holds. Since the Brownian motion has a.s. a 5-modulus of continuity 
bounded by 5 K for any k < 1/2, we choose k — 1/4 and take e proportional to 5 l / A . 
The maximal length M of the chains is fixed, large enough so that M>1 + ^ = 17 (see 
(jHJ)). Taking a unique scale parameter m G N we thus choose 

^(m) = (l + 3r + )M2 4 "\ 5{m) = ^. (4) 

It will be clear in ( ITT]) why this, with a suitable e(m) = pf, is a right choice. 
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We now define 

= :X« y(w )^(^^)n^(J-,M|)}. (5) 

We show in the next proposition that the set Q y is of full measure with respect to any hard 
globule Poisson process. 

Proposition 3.6 For any hard globule Poisson process \x G IT 9; one has 

[ P(Q y ) n(dy) = 1. 
Jm 

As a corollary, for \i a.e. y,P(£l y ) = 1. 

Proof We have to prove that / P(Q y ) ^{dy) = 0. 

JA g 

Thanks to Borel-Cantelli lemma, f A P(Q y ) fJ>(dy) vanishes as soon as the series 

and £ m J A P { 3k ^ : X ' (m) * y (^r) e ^(^M^) 

converges. Since for large £,Q e,y and the law of X ' y with initial distribution i^(-\yB(o,e) c ) 
are close, a similar argument as in [3] Proof of Proposition 3.2 and the condition (TTJ) yield 
that these series converge as soon as 

E m [ Q e{m ^(^:w((X ll X t ,^)>^)^dy)<+oc (6) 

J Ag 

and E m ^ ^ (m) ' y (3fc<2 4m :X(^)GGD M (^)) / ,(rfy)<+oo. (7) 
Following Proposition 13.21 

QW*(3i : w{{X h X u ^) > i) < Cl £(m) 3 2 4m exp (-c2 2m ) 

< c 4 2 16m exp (-c2 2m ) 

for a certain constant C4 > independent of y. The above right side is the general term of 
a summable series in m. Therefore (jSJ) holds. 
Following Proposition 13.51 




for a certain constant C5 > independent of y. We chose M large enough to ensure the 
summability in m of the right side, so (jJJ) holds and the proof is complete. ■ 
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3.3 The convergence 

In this subsection, y £ A g is still fixed and we study the convergence of the approximating 
processes as I — > +00. 

Proposition 3.7 For every uj in Q y and every i £ N, the sequence 

((X? m) * ^ (m) ^(u^ G N,t £ [0,l]) mGN , is sta- 

tionary as an element of C([0, 1], 1R 3 x R x x T/ie /z'mz't u>z7/ 6e denoted by 

((X?*,X?*)(u,t), L% y (u,t),L% y (oo,t),LrJ y (oj,t),j £ N,t £ [0,1]). Tfe/ore, 

lim X e(m) ' y (uj,-) = X°°' y (u,-) 

m—t+oo 

inC([0,l],Ag). 

Proof The main idea is that if a fixed globule moves along a nice path, it will only collide 
into a finite number of other globules. Thus dynamics (£ g ) reduces to an infinite number 
of SDE involving only a finite random number of particles up to time 1. 

Take u £ fi y and p > 0. Then, for m large enough, X^ m ^ y (u) and X^ m+1 ^ y (uj) both 

1 2 3 =1 2 6 
belong to the same set of regular paths J\f(—^—, — ) D Mi—r- , M, — ). 

For X in this set and k — 0, • • ■ , 2 4m — 1, we define the finite set of indices Jfc. m (X) as: 
J fc , m (X) := |i £ N, )| < v k>m or 

Xi{-^^) belongs to some — — chain of globules which intersects B(0, ffc jm )|, 
where v k , m := p + (1 + 3r + M)2 4m - 3r + Mk. 

One can show (similarly as in Lemma 3.3 [5]) that, for m large enough: 

{z : \Xi(0)\ <p}c J 2 ^_ 1)m (X) C • ■ ■ C J ,m(X) (9) 
Moreover, a globule with index in J^ m {X) does not bump into globules outside this set : 

• £ J fcim (X), j ? J fc , m (X) =► Vt £ [A ^±i] - > + + (10) 

and it stays in a large ball around the origin : 

z £ J fc , m (X) Vt £ [^L, f£±] |X,(t)| < z; fc _ lim < £(m) - 2r+. (11) 

Thus the penalization functions ^ m ^ y and ^ m+1 ^ y vanish on globule (X^t) , X(t)) if 
i £ J fc , m (X) and t £ [5^,5^]. Consequently, the paths X e ^' y (u) and X^ m+1 )> y (u;) 
satisfy the following simplified version of equation (£^' y ) : 

VA;£{0,---,2 4m -l} Vz£j fe , m (X) Vtef^,^], 

X i (t) = X i (-^-)+Wi(w,t)-W i (a;,-^-)+ V f 5^dL«(s) 

jGJ fe , m (X) 2^ n y 

x,(t) = x{ A) + t) _ ff wi( Wj A.) 

£ (^(t) - ^(^)) - (Li+(t) - M^)) + (M^ - ^-(^)) ■ 
ie^,m(x) 
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The initial configurations X^ m )' y (u;, 0) and X £< - m+1 )' y (w, 0) are equal to the same configu- 
ration y. Hence the set of indices J 0)m (X £ ^' y (u;)) and J 0im (X^ m+1 )' y (u;)) are equal and 
(xf m),y (w),xf m),y (w),z G J ,m.(X^ m )' y (u;)) satisfy the same equation as 
(Xf m+1) ' y (a;),xf m+1) ' y (u;),z G J , m (X^ m+1 )' y (u;)) during the time interval [0; The 
strong uniqueness in Proposition 13. II implies the equality of the final values X.^ m ^' y (u, ^r) 

and X^ +1 " y (a;, ^r) for the indices i in the set J 0>m (X^ m ) ,y (c<;)), which contains both sets 
Ji. m (X £ ( m ) ,y (u;)) and Ji im (X^ m+1 ) ,y (a;)). Thus these two sets of indices are equal, which 
in turn implies that the paths {xf m) ' y (u), Xf { m) ' y (u)) and (X? m+1) ' y (u;), X? m+1) ' y (u;)) 
coincide up to time for indices i G Ji jm (X^' m ) ,y (cu)) . Using inclusions ((HJ and the 
strong uniqueness again, we obtain the equality of both paths on J 2im (X^ m ) ,y (to)) = 
J 2)m (X £ ( m+1 )' y (cj)) up to time ^r, and so on. 

Strong uniqueness of the solution of (£^' y ) holds for the path X and the reflection term 
(linear combination of local times), but a priori not for each local time separately. How- 
ever, as shown in the proof of corollary 3.6 in [3J, the local times Ly, L i+ , Lj can be chosen 
in a unique way. With this choice, the same argument as above prove that local times 
L^iu^^t^ico^Lt^i^t) and L£ m+1) '>, t), L^ y (u, t), L? +1) ' y (u;, t) co- 
incide for i,j in J 0im (X £ ^' y (a;)) and t G [0; ^r], and then again for i,j in Ji )Tn ,(X^ m -'' y (a;)) 
and t < ^r, and so on. 

In particular, if \yi\ < p, then for m large enough depending on p, X? m) '» = xf m+1) ' y (u;) 
and xf m) '>) = If m+1) ' y (u) on the whole time interval [0; 1]. The associated local times 
can be chosen in such a way that they also coincide, which implies the equality of the 
reflection terms up to time 1. This completes the proof of the stationarity of the sequence 
of continuous functions ((Xp m '' y , Xf- m ^' y , Lfj 1 ' ,y , y , Lf^' y )(u, -))meN* and therefore 

its convergence to some path denoted by (X°°' y , X°°' y , L°°' y , L°^ y , L°^' y )(w, ■). 

To check the convergence of X £ ^ ,y (u, ■) in C([0, 1], A g ), we remark that for each continuous 

function / on M 3 x R with compact support, 

< X^)< y (u;, ■),/>=£ f( X ? m) ' y (^ ^ M '>, •)) 

i 

where the sum is indeed finite due to the minimal distance 2r_ between any pair of points 
Xf- m ^ ,y (oJ,t) and X^~ m ^ y (uj,t) and the local boundedness of path oscillations. Therefore 
the stationary convergence of each term insures the stationary convergence of the sum. ■ 

3.4 Properties of the limit process 

To complete the proof of Theorem 12.21 it suffices to show the following proposition. 

Proposition 3.8 For every y G A g := {x G A g : -P(f2 x ) — 1} the family of processes 

(Xp y (t),Xp y (t),L% y (t),L% y (fyL°^ y (t),i,j G N,t G [0,1]) with initial configuration y 
solves uniquely the stochastic equation (S g ). 

Proof The equation satisfied by each X°°' y (resp. X°°' y ) includes by construction a finite 
sum of local time terms. It is straightforward to prove, using a similar argumentation as 
in Proposition 4.1 of [5], that in fact this finite sum is already equal to the infinite sum 
present in (S g ). 
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The stationary convergence in Proposition 13.71 implies that: 

V/CN finite , 3m , Vm > m (X°°' y , X?*)^ G , ^) H A^, M, ^) 

Consequently, using the strong uniqueness in Proposition 13.11 as in the proof of Proposition 
13.71 uniqueness of the solution of (E g ) can be proved in the path space 

{XeC([0,lU 9 ): V/ CN finite, Vm , 3m >m (X h G Af^, i-)oA/(^L M, | 



Let us conclude with the proof of Proposition 12.31 that is with the reversibility of the 
solution of (£ g ) when the initial distribution is a hard globule Poisson process. 



Proof of Proposition 12.31 Using similar estimates as in the proof of Proposition 13.61 the 
solution of {Eg) starting with a hard globule Poisson process is the limit of processes whose 
distribution are close to Q e,y , which is a time-reversible measure. More precisely, we have to 
prove that, if \i G for any fi, . . . , fk bounded continuous functions on A g with compact 
support and for t\, . . . , i& G [0, 1] 

k k 

/ / * IfiiX^yfaU)) P(du) fi(dy) = ' iMX^yfal-U)) P(du) fi(dy) (12) 

which is equivalent to 

/ (n/i( x ' (m)jr (*i))-II^ x ' (m),y ( 1 -**))) dp ^ d y) = Q - 

- l s J n i=l i=l 

By computations similar to those done in [5] to obtain inequality (17), we have 



k k 

I f (n fip&b)) - n ^( x '' y (! - **))) dp M^y) 

J A g Jn i=1 i=1 

k k 

H fi(X(U)) - II - **))) Q e ' y (dX) fi(dy) 

i=i 

+2 J] sup |/ 4 (x) 



A •'ft i=l 



i=1 XGA 



1 — ) fx(dy) 



The first term of the right hand side is equal to 0. The second term tends to zero as I tends 
to infinity, thanks to assumption ([!]). ■ 
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